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Abstract 

We review what is known about the contributions to the final entropy from the diff'erent stages 
of a relativistic nuclear collision, including recent results on the decoherence entropy and the 
entropy produced during the hydrodynamic phase by viscous efl'ects. We then present a general 
framework, based on the Husimi distribution function, for the calculation of entropy growth in 
quantum field theories, which is applicable to the earliest ("glasma") phase of the collision during 
which most of the entropy is generated. The entropy calculated from the Husimi distribution ex- 
hibits linear growth when the quantum field contains unstable modes and is asymptotically equal 
to the Kolmogorov-Sinai (KS) entropy. We outline how the approach can be used to investigate 
the problem of entropy production in a relativistic heavy-ion reaction from first principles. 



1. Overview 

The agreement of hydrodynamic calculations of the flow anisotropy of the matter produced 
in nuclear collisions at RHIC with the elliptic flow measurements rests on the assumption of an 
early equilibration of the matter on a time-scale of the order of 1 fm/c. An important problem in 
the description of relativistic heavy ion reactions is thus to understand how the produced matter 
equilibrates so quickly. From a thermodynamic standpoint, this question can be answered by 
studying when and how entropy is created in the reaction. One can distinguish five different 
stages of entropy production: 

1. Decoherence of the initial nuclear wave functions; 

2. Thermalization of the partonic plasma ("glasma"); 

3. Dissipation due to shear viscosity in the hydrodynamic expansion; 

4. Hadronization accompanied by large bulk viscosity; 

5. Viscous hadronic freeze-out. 

What is known about the contribution of the different stages to the final entropy is depicted in 
Fig. [T] First, we note that the final entropy per unit rapidity dS/dy is one of the best known 
quantities in relativistic heavy ion physics. In the case of Au-i-Au collisions at RHIC dS /dy at 
freeze-out can be determined from an analysis of the final hadron spectra in combination with 
the information on the source radius derived from identical particle (HBT) correlations 1 1 1. The 
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slightly extrapolated result for the 6% most central Au+Au collisions at -\fs^ - 200 GeV is 
(dS /dy)f - 5600 + 500 at midrapidity |2j. Alternatively, the final entropy can be deduced from 
the measured hadron abundances, combined with the calculated entropy per particle for a hadron 
gas in chemical equilibrium, which yields the result (dS /dy)ch = 5100 + 200 |^. The 10 percent 
diff'erence can be attributed to the entropy production during the hadronic freeze-out and reflects 
the significant viscosity of a hadronic gas 131. 
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Figure 1 : Entropic history of a central Au+Au collision at top RHIC energy. The values of dS /dy indicate the entropy 
per unit rapidity reached at the end of various collision stages based on experimental data and model estimates. 



2. Mechanisms of entropy production 

Broadly speaking, we can distinguish between entropy production by decoherence of the 
initial quantum state, by the equilibration process, and by dissipative processes in the hydrody- 
namic expansion phase. Starting with decoherence, the loss of coherence is measured by the 
decay in off-diagonal elements of the density matrix p. A practical way of investigating this is 
to calculate the decay rate of the quantity Trp^/[Trp]^. The time scale of the decoherence of the 
initial nuclear wavefunction has been studied by this method in the color glass condensate model 
with the result that the characteristic decoherence time is |j4j|5l Tdeco = cQ^', where Qs is the 
nuclear saturation scale, and c denotes a calculable constant close to unity. While this result is 
expected on dimensional grounds, the fact that c < 1 is important because it shows that deco- 
herence occurs over a time of less than 0.2 fm/c. The amount of entropy created from loss of 
coherence alone can be estimated by multiplying the number of causally disconnected transverse 
domains in a central Au-i-Au collision by the entropy Sdeco ~ (inlnh + l)/2 obtained from the 
complete decoherence of a coherent quantum state with average occupancy «. Accounting for a 
longitudinal coherence length Ay ^ I /as and using h x I /a^ one obtains Ii5ji6jj: 

(dS/dy)d,co * ^e?^'ff.(ln27r/a, + 1) * 1,500 (1) 
2 



at midrapidity for a central Au+Au collision, roughly one quarter of the final entropy. 

The other stage which is reasonably well understood is the hydrodynamic expansion stage. 
RHIC data on the elliptic flow of hadrons in non-central collisions indicate that the shear viscosity 
during this phase is small. The bounds on the shear viscosity rj compatible with the RHIC data 
lie characteristically in the range 77/s < (2 - J)! An ll7ll8ll9l [T0l[Tll . Such a small viscosity cannot 
contribute much to the entropy of the expanding fluid. If one sets the start of the hydrodynamic 
era at Tth ~ 1 fm/c, where Tjh denotes the thermalization time, one-dimensional solutions of boost 
invariant hydrodynamics show that the the entropy of the fluid increases by about 15 percent 
between Tth and the time of hadronization lfT2l . This implies that the entropy at thermalization is 
approximately (dS/dy)th ~ 4,500. 

3. Entropy growth rate 

These considerations tell us that at least half of the final entropy must be generated during the 
thermalization process. If the quark-gluon plasma were a weakly coupled system, the growth of 
the entropy during this stage could be calculated from a partonic Boltzmann equation. However, 
the small shear entropy and other observations from the RHIC experiments, such as jet quench- 
ing, indicate that the plasma is strongly coupled. In addition, it is thought that the initial state (the 
"glasma" 1 13 , 14 1), as well as the pre-equilibrium quark-gluon plasma fTSl, are characterized by 
strong, partially coherent color fields. Any formalism that hopes to describe entropy production 
during this early stage must, therefore, be able to describe the growing complexity of a quantum 
system making the transition from a regime of field dominance to the hydrodynamic regime. We 
have recently proposed such a general formalism [16 J . 

The idea is basically simple: the growing entropy measures the increasing intrinsic complex- 
ity of the quantum state of the system after appropriate coarse graining. The problem is how 
to impose a minimal amount of course graining without assuming the answer The solution to 
this problem dates back to Husimi IfTTl . The Husimi distribution is a convolution of the Wigner 
function W{p, x, t) of the system with a minimum-uncertainty Gaussian wave packet: 



Here x and p stand for all "position" and " momentum" variables characterizing the system, 
which may include particle positions as well as field amplitudes (see ifTSl for the definition of 
the Wigner functional for quantum fields). The Husimi distribution depends on the squeezing 
parameter A, which measures the ratio of position and momentum uncertainty. Its value can 
be chosen at liberty, but once fixed, will not change with time. The Husimi distribution can be 
understood as a coarse grained phase space distribution of the quantum system, where the coarse 
graining corresponds to the projection on a coherent state of the quantum system. We recall 
that coherent states are the closest quantum analogues of classical systems compatible with the 
uncertainty relation. Because the Husimi distribution can be shown to be positive (semi-)definite, 
it permits the definition of a coarse-grained entropy, first introduced by Wehrl I.19J : 



Quantum systems containing unstable modes, i.e. modes with an exponentially growing ampU- 
tude, have a linearly growing Husimi -Wehrl entropy |161. The entropy growth rate is indepen- 
dent of the squeezing parameter A and given by the sum of the exponential growth rates of aU 
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unstable modes. In classical dynamical systems, this quantity is known as the Kolmogorov-Sinai 
entropy, or KS entropy, and defined as the sum over all positive Lyapunov exponents A/^ of the 
system: dSn.^ldt — > ^ks = Yik '^k 9{M)- The KS entropy is understood to be a measure of the 
growth rate of the coarse grained entropy of a dynamical system starting from a configuration far 
away from equilibrium, after an initial start-up phase during which unstable fluctuations grow to 
dominance and before it gets too close to its micro-canonical equilibrium 1201 . 

4. Entropy growth in QCD 

The Husimi-Wehrl entropy provides the basis for the formulation of a comprehensive ap- 
proach to entropy production in heavy-ion collisions. The lattice regularized Yang-Mills equa- 
tions describing the dynamics of classical color fields are known to be strongly chaotic 112111221 . 
and the KS-entropy of the classical Yang-Mills field was shown to be a thermodynamically ex- 
tensive quantity [23 1. The lattice approach can be extended to include Gaussian fluctuations 
ll24l . A calculation of entropy production using this approach with a random initial conditions 
is presently underway ^251, and it would be interesting to evaluate the evolution with the initial 
quantum fluctuations around classical glue fields in the colliding nuclei ||26l . 
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